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A polarized hadron collider experiment must have a great discovery potential for a 
search of physics beyond the standard model. Experimental data of various symmetry tests 
at RHIC are going to be obtained within a few years. The author developed a simulation 
tool, studying a sensitivity of hunting contact interaction at RHIC by measuring parity 
violating spin asymmetries. 



1 Introduction 

One of the most vital tasks of an experiment using a polarized hadron collider 
in a new energy region is the study of physics beyond the standard model. In the 
past few years, several theoretical works |[ have been devoted to the studies 
of the contact interaction (or, compositeness) M and new gauge bosons W'± § 
and Z' In these theoretical works, B it has been shown that RHIC can reach 
a similar sensitivity to that of the TEVATRON, due to its polarized beam. 

The purpose of the present study is to explore the discovery potential for physics 
beyond the standard model from the experimental perspective. Considering the 
family pattern of quarks and leptons, compositeness is an fascinating idea. For 
there are no standard sub-quark models, model independent treatment using a 
compositeness scale A is often used. We call such new generalized interaction as 
the contact interaction. One possible Feynmann diagram is shown in Fig.|l|. For 
any kind of an new interaction with energy scale A can be taken into account, for 
example, new gauge boson exchange between "standard" quarks can be understood 
as a kind of the contact interaction. A is a model-independent scale parameter of 
the contact interaction, defined as F(Q 2 ) = (1 + Q 2 /A 2 )~ 1 , where F(Q 2 ) is a 'form 
factor' of the quark and lepton. Therefore, sometimes one interpret the physical 
meanings of the the contact interaction as quark (lepton) compositeness. However, 
the contact interaction has general formalism, so any kind of new interaction can 
be roughly included. 



2 Simulation Procedure using PYTHIA+POLBY 

In order to study physics sensitivity and to determine a search window, sim- 
ulation studies with an event generator are indispensable. Several non-standard 
scenarios can be examined using PYTHIA. The contact interaction can be also ex- 
amined through it. However, as with all other subprocesses, PYTHIA includes only 
helicity-averaged cross sections for the contact interaction. The author developed 
a plug-in program 'POLBY' for PYTHIA by which hadronic spin asymmetries for 
Drell-Yan and quark-scattering process can be examined. Helicity-selected matrix 
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elements have been theoretically obtained in LO not only for the standard model 
but also for the contact interaction. In case of Drell-Yan process, corresponding 
matrix elements were obtained from Virey's calculation for DIS by crossing A 
weighted mean method is applied in order to extract spin asymmetries using "unpo- 
larized" event generator PYTHIA. A weight factor is calculated for each generated 
event, then its mean value of the weight factor distribution gives us the final spin 
asymmetries. First, partonic cross section a can be obtained at a kinematics region 
by summing up the matrix elements for all the possible interactions. In order to 
estimate the hadronic spin asymmetries, we need hadronic cross sections a hlfl2 with 
helicity hi, hi selection. They can be obtained by combining with polarized quark 
distribution functions . For example, cr ++ can be obtained as; 

cr ++ = <7 ++ <7+<7+ + 0- - ~ft~?2 + a+-q+q^ + a^+q^q}. (1) 

Using the hadronic cross sections, we can estimate the weight factor in combinations 
between different spin selected hadronic cross sections. For single spin asymmetry 

er~ — <j + a + a — a ++ — 
L = a- + cr+ = a— + <j-+ + a++ + a+~ () 

and parity violating double spin asymmetry 

All - a — ^ 0) 

corresponding weight factors can be obtained by making the same cross section 
ratio in event-by-event. The distribution of the weight factor depends on the quark 
distribution functions, therefore, the process of making a mean value on these dis- 
tribution implies an Monte-Carlo integration over the quark distribution functions. 
There are two free parameters in the matrix elements for the contact interaction. 
First one is the scale parameter A and the second one is a sign of the interference 
between the standard model and the contact interaction. 
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3 Drell-Yan process 

Parity violating spin asymmetries are examined using PYTHIA+POLBY. In 
PYTHIA, Drell-Yan events can be generated in single subprocess, which include 
the standard model 7*/Z° exchange and also the contact interaction with its scale 
A. Fig. U shows the results on Drell-Yan yield with selections of beam energies (-y/s 
= 500, 750, 1000 GeV) and the scale A. Here, GS-A polarized parton distribution 
function || is used. At any beam energy, excess from the standard model are 
observed when we require A = 1 TeV, however, the excess is negligible if A > 2 
TeV. Considering the percent-level systematic error in a cross section measurement, 
observation of the cross section excess is hopeless. In Drell-Yan process, contribution 
from the pure contact interaction is very small because of its A dependence of A~ 4 , 
on the other hand, the interference terms between the contact interaction and Z°, 
7* have A -2 dependence. 

On the other hand, a polarized experiment has a great advantage because of 
its systematic error cancellation. Simulation results using PYTHIA+POLBY on the 
parity violating single spin asymmetry Al are shown in Fig. |§|. Pure Z° exchange 
and 7* — Z° interference contribute on the parity violation on the standard model. 
We can sec the large deviation from the standard model if we include the contact 
interaction. Dominant contribution on the parity violation is the interference term 
between the contact interaction and the standard model. In FigJ^, we can see a clear 
deviation from the standard model even at A = 5 TeV, however, the experimental 
statistical error have to be also large. Because of its relatively small cross section, 
required integrated luminosity is large as shown in Fig.^J. Considering the RHIC 
luminosity of 800 pb^ 1 at y/s — 500GeV, our sensitivity is limited around A ~ 1 
TeV, which is lower than the already established limits in e + e~ collisions for the 
electron-quark contact interaction (A > 5.4 TeV) by ALEPH and in pp collisions 
for the muon-quark contact interaction (A > 2.9 TeV) by CDF JlO[ j . 

4 Inclusive Jet production process 

Although Drell-Yan process is simple and clean, we need to examine other pro- 
cess, which has large production cross section because of the small yield on Drell- 
Yan process. Parity violation in the inclusive jet production process has already 
examined by J.-M. Virey theoretically [Q. The jet production cross section is shown 
in Fig.|5| as functions of jet Pt- Lower Pt region is dominated by gluon-related pro- 
cesses, on the other hand, higher Pt region is dominated by quark-quark scattering. 
The contact interaction couples to leptons and quarks, but not to gluons in the cur- 
rent model. It means that we do not think gluon is a composite particle. Therefore, 
sensitive Pt region of the contact interaction must be limited to the higher Pt 
region. In case of Drell-Yan process, all the interactions we need to consider, are 
the standard model 7* /Z°, the contact interaction, and their interferences. On the 
other hand, we need to take into consideration QCD, 7*, Z°, W^, the contact in- 
teraction, and their complex interference terms for the jet production. In PYTHIA, 
above processes must be generated as independent subprocesses. This restriction 
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cause severe problem when we consider the interference terms between different 
subprocesses in PYTHIA. It is because there cannot exist interference effects if we 
distinguish the subprocess. In case of Drell-Yan process, we can use single subpro- 
cess, therefore, the interference effects can be automatically included. In order to 
restore quantum mechanics in the jet production, a special treatment is required. 
In POLBY, a hadronic spin asymmetry A is estimated using event- by- event asym- 
metry weight factor W / ( J 4); 

^^i— event v^) /^\ 

^-^i— event 

Here, the "events" are generated by PYTHIA subprocess. The "interference correc- 
tion" can be done by 

^process 

W* (A) - W* (A) X ^ocess=QCD,Weak EM,CTV (g) 

with Eq.|]. Then if we use pure QCD subprocess for the PYTHIA event generation, 
contribution from other processes on the asymmetry can be taken into account 
by summing up over all the processes. Fig.^J shows a sample distribution of W(A) 
before/after the interference correction of EqJ|. The mean value of the IF (A) dis- 
tribution shows the final hadronic spin asymmetry. By the operation Eq.||, W (A) 
distribution is distorted and we can get the corrected asymmetry from the distorted 
distribution. By using Eq.|^, we can estimate accurate asymmetries. For an exam- 
ple, parity violation around Jacobean peak is examined. In Fig.|?], we can confirm 
the large contribution from the interference term on the parity violation around 
the Jacobean peak. Subprocess-separated plot in Fig.f7| shows that qq' — ► qq' sub- 
process dominates the parity violation at the W peak. Fig|| shows a preliminary 
results, which can be compared with the theoretical calculation Q|. Considering 
the difference on the pseudo rapidity cut and binding, the results meets well with 
the calculation by Virey. Using these results, discovery potential at RHIC can be 
explored. In Figjp, required integrated luminosities are plotted as functions of A. 
Here combined results from the present study at lower A region and that from Virey 
at higher A region are shown. It is shown that we can reach more than 3 TeV 
sensitivity, which is higher than the current limit of 2.7 TeV reported by DO for the 
quark-quark contact interaction, jlo| with the current RHIC plan. It is also shown 
that the sensitivity limit can be more than 6 TeV using the RHIC upgrade plan. 
We can calculate many combinations of parity violating spin asymmetries using 
PYTHIA+POLBY with many realistic experimental conditions. For the procedure 
using weighted mean can be understood as a kind of Monte-Carlo integration, the 
calculation speed is very fast. Therefore this is a strong tool to examine an ex- 
perimental sensitivity and also to make a standard model reference after getting 
experimental data. 

However, in order to use it as a standard model reference, NLO calculation is 
indispensable. NLO effect is difficult to take into an event generator. The situation 
is similar to the treatment of the interference. One large assumption for the inter- 
ference correction is that, pure QCD process (generated events) and other processes 
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must have their same final states. In principle, NLO correction using similar cross 
section correction like; 

^process 

W* (A) -> W* (A) X ^Process=NLO+LO (g) 

<J 

should be able to be applied. Theoretical calculations in NLO are now under 
progress We should and will be able to be ready to examine experimental 
data from RHIC in near future. 
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Fig. 1. Possible Feynmann Diagram of the contact interaction with energy scale A. 
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Fig. 2. Cross section of Drell-Yan process is shown including the contact interaction. 
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Fig. 3. Single spin asymmetries At, are drawn as functions of Drell-Yan mass with selec- 
tions of A. 
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Fig. 4. Required integrated luminosity as functions of A is shown. 
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Fig. 5. Jet production cross section is shown. At high Pt region, quark-quark scattering 
is dominant. Small bump around Mw/2 corresponds Jacobean Peak 




Fig. 6. Sample plot of the weight factor. Hatched(Blank) histogram shows before(after) 

the interference correction. 
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Fig. 7. Fine structure of the Jacobean peak on parity violating spin asymmetries. 
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Fig. 8. Expected parity violating spin asymmetries. 
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Fig. 9. Required integrated luminosities are plotted as functions of A. 
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